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Abstract 

We introduce suitable coordinate systems for interacting many-body systems 
with invariant manifolds. These are Cartesian in coordinate and momentum 
space and chosen such that several components are identically zero for motion 
on the invariant manifold. In this sense these coordinates are collective. We 
make a connection to Zickendraht's collective coordinates and present certain 
configurations of few-body systems where rotations and vibrations decouple 
from single-particle motion. These configurations do not depend on details of 
the interaction. 
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I. INTRODUCTION 



Dynamical systems with invariant manifolds in phase space have been the subject of 
ongoing research in recent years. Many authors have considered the case of two or more 
coupled identical systems that are chaotic. On invariant manifolds the subsystems display 
identical or synchronized motion, and the manifold's transverse stability is investigated |T] [7J . 

An alternative approach is based on the observation that any rotationally invariant sys- 
tem of identical interacting particles possesses low dimensional invariant manifolds in clas- 
sical phase space ||. On such manifolds, the classical motion displays largely collective 
behavior and decouples from more complex single-particle behavior. The importance of a 
given invariant manifold depends crucially on its stability properties. If the manifold under 
consideration is sufficiently stable in transverse directions, the quantum system may exhibit 
wave function scarring PHTll or display a strong revival for wave-packets localized to the 



vicinity of the manifold [12] . These findings may be directly associated with the slow decay 
of collective motion due to the coupling between collective and single-particle motion. 

In this paper we propose suitably adapted coordinate systems that separate collective 
and single-particle motion on the invariant manifolds mentioned above. Such coordinates 
clarify the separation of collective and single particle motion and may be useful in several 
applications. We have in mind (i) the problem of damping and dissipation of collective 
excitations and the interplay of collective and chaotic motion in atomic nuclei [|T^ -pO|, which 
is often addressed in the framework of single-particle motion in a time dependent mean field; 
(ii) multi particle fragmentation of atoms at threshold which evolves over highly symmetric 
configurations corresponding to invariant manifolds @j2l|]; (iii) the structural stability of 
invariant manifolds |22| . 



There is a traditional way to introduce collective and single-particle coordinates in inter- 
acting many-body systems. Aiming at the description of nuclear vibrations and rotations, 
Zickendraht |23[ introduced a system of collective coordinates in a self-bound many-body 
system. Three of these coordinates describe the center of mass motion, and six collective 
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degrees of freedom govern the dynamics of the inertia ellipsoid. The remaining coordinates 
are of single-particle nature. We shall establish the relation of coordinates of invariant man- 
ifolds to those defined by Zickendraht. Furthermore we shall show that more complicated 
collective motion, e.g. shearing modes can be described 

This article is divided as follows. In the next section we introduce suitable coordinate 
systems for interacting many-body systems with invariant manifolds. We give a construction 
recipe and present a detailed example calculation. As an application we give a potential 
expansion around an invariant manifold and discuss stability properties. In section |IV| we 
make a connection with the Zickendraht coordinates. We present examples where the motion 
of the inertia ellipsoid corresponds to the motion on an invariant manifold. For such initial 
conditions the traditional collective motion decouples completely from the single-particle 
degrees of freedom. We also find that collective coordinates as defined here are capable of 
other types of motion. Therefore we finally discuss how motion on or near such invariant 
manifolds could be interpreted as collective motion of a system. 

II. COORDINATES FOR INVARIANT MANIFOLDS 

In this section we present a transformation from Cartesian single particle coordinates 
in position and momentum space to Cartesian coordinates that are adapted to invariant 
manifolds. The new coordinates consist of "collective coordinates" that govern the motion 
on the invariant manifold and of coordinates transversal to this manifold that represent the 
single-particle aspects. 

Consider rotationally invariant systems of N identical particles in d spatial dimensions 
(d — 2 or d — 3). The Hamiltonian is invariant under both, the action of the rotation group 
0(d) and the group of permutations Sjv- One may now take a finite subgroup Q C 0(d) 
with elements g and properly chosen permutations P(g) such that 

9P' 1 (9)(P,Q) = (f,q), Vg^G p = (p u ...,p Nd ),q=(q 1} ...,q Nd ) (1) 

for points (p, q) on some invariant submanifold of phase space. On such a manifold, the 
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action of certain rotations g can be canceled by permutations. These permutations clearly 
form a subgroup isomorphic to Q. 

Fig. [1] shows a configuration of four particles in two spatial dimensions that corresponds 
to a point on an invariant manifold. The operations of elements from the discrete symmetry 
group Q = C 2v can be undone by suitable permutations of particles. This leads to a collective 
motion with two degrees of freedom which we shall identify with vibrations. 

Fig. |] shows two spatial configurations of eight (pja) and six @y) particles, respectively 
which display a symmetry. If initial momenta display the same symmetry the motion on 
the invariant manifold will have two degrees of freedom. For eight particles the radii of the 
two circles will oscillate synchronously, and the two circles will vibrate against each other. 
For the six particles we will have a vibration of the radius of the circle and of the two particles 
along the vertical axis. We may choose initial momenta to reduce the symmetry group to 
C&h which will allow rotations around the vertical axis and thus add an additional degree of 
freedom. For eight particles we could alternatively choose initial conditions that are limited 
to a -D4 symmetry. Besides the vibrations discussed above this would allow for a shearing 
motion of the two circles thus yielding again three degrees of freedom. We could also reduce 
the fourfold rotation axis to a twofold one and obtain D 2 h, D 2 or C 2v as remaining symmetry 
groups yielding more collective degrees of freedom. Adding two particles symmetrically onto 
the principal axis of rotation would also increase the number of degrees of freedom by one. 
Other reductions of symmetry will yield different invariant manifolds with varying degrees 
of freedom. We will see this exemplified by explicit construction of coordinates. 

We may use the definition ([I]) directly for the construction of coordinate systems where 
invariant manifolds correspond to coordinate axis or planes, i.e. non collective coordinates 
vanish for motion on the invariant manifold. To this purpose we consider the many-body 
system in Cartesian coordinates in momentum and position space. In what follows we will 
introduce orthogonal transformations in configuration space only; momenta will be subject 
to the same transformation. 

In a Cartesian coordinate system each element g G Q and each permutation P(g) can 
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be represented by an orthogonal matrix M 9 and P 9 of dimension Nd. It is clear that 
the products M^Pj form a matrix group TC that acts onto position and momentum space, 
respectively. The construction of the coordinate system is now straightforward. Every vector 
p and q may be expanded in basis vectors of the irreducible representations (IRs) of 7i by 
means of projectors . 



n„ = ExS"HPj- (2) 

geg 

Here Xg denotes the character of g in the z/th IR. Similar formulae hold for momentum 
space. The projection onto the identical IR defines the invariant manifold. Note that the 
identical representation is one-dimensional while the invariant manifolds of interest typically 
have higher dimensionality. We can find independent vectors on the manifold by projecting 
from different vectors, but in practice the construction of the independent vectors seems to 
be unproblematic as we shall see in the example. 

A comment on the rotation symmetry is in order. Like any Cartesian coordinates, the 
coordinates introduced in this article do not explicitly reflect the invariance under rotations. 
Acting on an invariant manifold, rotations generate a continuous family of equivalent man- 
ifolds. Our coordinates, however, single out one particular manifold. For quantum systems, 
the rotation operator may easily be constructed and used for projection onto subspaces of 
definite angular momentum. 



III. A SIMPLE EXAMPLE 

We now illustrate the proposed construction explicitly for four particles in two dimen- 
sions and a quartic potential, considering the invariant manifold shown in Fig. |l[ We shall 
also expand the potential near the invariant manifold to second order in the transversal 
coordinates. 

The invariant manifold is defined by those points which are invariant under 7i = 
{E, u x P(i2)(34), Ci/-P(i4)(23), C2-P(i3)(24) }, where E denotes the identity, P a permutation of 
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particles as indicated, a a reflection at the axis indicated, and C 2 a rotation about ir. Thus, 



H = C 2v with four IRs labeled by v = A u B x , A 2 , B 2 @. Let q = (x 1} x 2 , x 3 , x 4 , j/i, y 2 , j/ 3 , j/ 4 ) 
denote a coordinate vector in position space (xi, yi denote the coordinates of the z'th parti- 
cle). We have 

Eq= (X!,X2,X3,X4, J/l, J/ 2 ,J/3, J/4), 
<TxP(l2) (34) q = (x 2 ,Xi,X4,X 3 , -J/ 2 , -J/1, -J/4, -J/3) 
C 2 P(13)(24) q= (-X3, -X 4 , -Xi, -X 2 , -J/3, -J/4, -J/i, -J/ 2 ) 
0"^(14)(23) <f = (-^4, -^3, -X 2 , "^1, J/4, J/3, J/2, J/l)- 

Using the character table of |2l| and the projectors (|2|) one constructs the following 
basis vectors corresponding to the IR labeled by 



A,: 


e[ = 


(1, 1,-1, -1,0,0, 0,0) /2, 


e 2 = (0,0,0,0,l,-l,-l,l)/2, 


B x : 


e's = 


(l,l,l,l,0,0,0,0)/2, e' 4 


= (0,0,0,0,l,-l,l,-l)/2, 


A 2 : 


4 = 


(1, -1,-1, 1,0,0, 0,0) /2, 


e' 6 = (0,0,0,0,l,l,-l,-l)/2, 


B 2 : 


e' 7 = 


(1, -1,1, -1,0,0, 0,0) /2, 


e' 8 = (0, 0,0, 0,1,1, 1,1) /2. 



The vectors associated with the identical IR A\ span the two-dimensional invariant manifold 
and the vectors associated with the IRs B x , A 2 , B 2 span the transverse directions. 

We now present the orthogonal transformation that transforms the single particle co- 
ordinates q into the coordinates adapted to the invariant manifold. In our example x and 
j/-components do not mix and we have 
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To illustrate the example and to further demonstrate the usefulness of the newly intro- 
duced coordinate system we want to consider the the interacting four-body system with 
Hamiltonian 



h = E ((pI + ?&)/2 + 16(^ + y 2 ) 2 ) - E [0* - + (vi - yj) 2 } ■ (4) 

i=l i<j 

This Hamiltonian has been studied previously ||12|| . In particular, the stability of the in- 
variant manifold displayed in Fig. [I] has been studied by computing the full phase space 
monodromy matrix of several periodic orbits that are inside the invariant manifold. It was 
found that several orbits are linearly stable in transverse directions or possess rather small 
stability exponents. Qualitatively, this may also be understood by studying the Hamil- 
tonian (Hj) close to the invariant manifold. We therefore use the transformation (|3]) and 
expand the potential of Hamiltonian to second order in the transverse directions labeled 
by (e 1; . . . e 6 ) while keeping the full dependence of the coordinates (x, y) inside the invariant 
manifold. One obtains the quadratic form e^\e where 
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A diagonalization of V yields the eigenvalues 



Ax 2 + I2y 2 



A 12 = 10x 2 + Uy 2 ± 2^x 4 + 25y 4 + 5Ax 2 y 2 , 
A 3i4 = 14x 2 + 10y 2 ± 2Jy 4 + 25x 4 + 5Ax 2 y 2 , 



A 5 6 = 6( x 2 + y 2 ) ± 2 J9x 4 - 2x 2 y 2 + 9y 4 . 



All eigenvalues are non- negative and vanish at the origin (x = y = 0). Thus, instability may 
occur only in its vicinity. Though the expansion of a potential around an invariant manifold 
is no substitute for the computations of Lyapunov exponents or monodromy matrices, it is 
a first step when estimating stability properties of such manifolds. 
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IV. ZICKENDRAHT'S COORDINATES AND INVARIANT MANIFOLDS 



Almost thirty years ago Zickendraht P3J introduced a set of collective coordinates to de- 
scribe nuclear vibrations and rotations, as well as their coupling with single particle motion. 
We shall discuss to what extent these coordinates correspond to the ones we introduced in 
the previous sections. On one hand this will allow to identify certain vibrational modes of 
a many-body system with invariant manifolds. On the other hand we shall also see that 
our procedure proposes collective movements that are not of the type described easily in 
Zickendraht 's coordinates. 

Following Zickendraht [23[ we write the coordinates r*j of the i th particle in the center of 



mass system as 

ri = sn y{ + s i2 y 2 + s i3 y 3 , i = 1, . . . , N (5) 

where the yi span the inertia ellipsoid and Sik are non-collective coordinates which for sim- 
plicity we shall call single-particle coordinates. The newly introduced coordinates $ and Sij 
are not independent. The constraints are 

Vi-Vj =ViVAj, i,j = 1,2,3 

N 

$>y = 0, j = 1,2,3 

N 

s ij s ik = fijki 3i k = 1; 2, 3. 

i=0 

The first six equations ensure the orthogonality and normalization of the principal axis of the 
inertia ellipsoid whereas the next three equations fix the origin at the center of mass system. 
The last six equations are orthogonality relations of the single-particle coordinates. In the 
center of mass system, one may therefore characterize the iV-body system by its inertia 
ellipsoid (e.g. three Euler angles of the principle axis and three moments of inertia) and 
3iV — 9 single particle coordinates. The moments of inertia /; are related to the coordinates 

Vx by 

h = m(yl+yl), I 2 = m(yl + yl), h = m(yl + y 2 2 ), (6) 



where m denotes the mass of the particles. 

It is interesting to determine those configurations, where the motion of the many-body 
system may be described in terms of the collective coordinates yi only. While such motion 
would be restricted to some invariant manifold in phase space it would not obviously be one 
of those defined by eq. (|I]). We may however determine invariant manifolds (|I]) such that 
the motion on the manifold changes only the inertia ellipsoid of the system and hence may 
be described entirely by Zickendraht's collective coordinates y^. Two necessary conditions 
for this a situation are easily stated. First, the number of coordinates on such invariant 
manifold may not exceed six in the general case and three in the case of pure vibrations. 
Second, every motion on such an invariant manifold has to change the inertia ellipsoid of 
the many-body system. 

For simplicity let us start with the a system of four particles in two spatial dimensions 
and the invariant manifold displayed in Fig. [I], i-e. 



and the momenta are chosen by replacing x — > p x ,y — > p y . Computation of the moments 
of inertia yield the collective Zickendraht coordinates yi = 2x, y% = 2y. On the invariant 
manifold the remaining coordinates are given by Sn = Su = S21 = —S22 = — S31 = — S32 = 
— S41 = S42 = 1/2. This shows that every motion on the invariant manifold only changes the 
moments of inertia and therefore decouples from the single-particle motion. 

We next consider the example of an eight-body system in three dimensions. Let 
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(7) 



denote a configuration restricted to the invariant manifold displayed in Fig. [2] (a) with 
symmetry. (The momenta are chosen by replacing x — > p x ,y — » p y , z — *> p z in eq.([?p.) The 
moments of inertia are l\ — l-i — 4m(s 2 + y 2 ) + 8mz 2 , 7 3 = 8m(x 2 + y 2 ) and yield collective 
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coordinates @ y\ = y\ = 4(x 2 + y 2 ),y 2 = 8z 2 . Since the inertia ellipsoid is symmetric we 
have a freedom in choosing two of its principle axis. Using 
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one obtains constant single-particles coordinates Sn = — s 3 i = s 5 i = —s n = s 2 2 = — s 42 = 

^62 — — ^82 = 1/2? ^13 = S23 = S33 = S43 = — S53 = — Sq3 = — S73 = — S 83 = 1/V8 for 

the motion on the invariant manifold. Thus, the single-particle motion decouples from 
the collective motion on the invariant manifold. Similar results hold for the six particle 
configuration displayed in Fig. 

It is also instructive to consider one counterexample. The configuration 



r 5 



x 

y 

z 

X 

-y 





-y 




—x 




y 


, r 2 = 


X 


, ^3 = 


-y 


, r 4 = 


—x 




z 




z 




z 




y 




—x 




-y 


, re = 


X 


, r 7 = 


y 


, r 8 = 


—x 




—z 




—z 




—z 



displays D 4 symmetry and differs from configuration (0) by a shearing motion. Like in the 
previous example, the moments of inertia are given by I\ — I2 — 4m(i 2 + y 2 ) + 8mz 2 , 13 = 
8m(x 2 + y 2 ) and the ellipsoid of inertia is symmetric. However, no choice of the principal 
axis allows to fulfill eqs. fl5|) with constant single-particle coordinates s^-. Therefore, single- 
particle degrees of freedom depend on collective degrees of freedom and a decoupling does 
not exist using Zickendraht's coordinate system. A decoupling is obtained by using the 
coordinates introduced in this work. However, the collective motion on the appropriate 
invariant manifold does not correspond to pure vibrations or rotations of the inertia ellipsoid. 
These findings are interesting e.g. in relation with with the magnetic dipole mode in nuclei 
26j] since this type of collective behavior is associated with a shearing motion. 
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V. DISCUSSION 



We constructed an orthogonal transformation that maps the Cartesian single particle 
coordinates of a many-body system to a new Cartesian coordinate system that distinguishes 
collective and single-particle motion. The collective degrees of freedom govern the motion 
that is restricted to a low- dimensional invariant manifold and are decoupled from single- 
particle degrees of freedom on this manifold. We have demonstrated that there are several 
configurations of few-body systems, where the motion on the invariant manifold corresponds 
to a vibration or rotation and may be described in terms of Zickendraht's collective coordi- 
nates, but differs when the collective motion goes beyond that. These results are indepen- 
dent of the details of the Hamiltonian of the iV-body system, and are entirely determined 
by rotational and permutational symmetry. 

Using the results of this article as well as those of refs. we can draw the following 

picture: First it is possible that an invariant manifold is spanned exactly by the vibrational 
and rotational modes of a few-body system; second such manifolds may be stable or have 
small instability exponents in transversal directions; third the revival probabilities of wave 
packets launched on such manifolds are large; last, as a conclusion of these points we may 
have a collective motion near the manifold whose damping is characterized by the decay rate 
in transversal direction. We also found that there may be other collective motions; this was 
displayed in an example of shearing motion, but there can be others such as breathing modes 
etc. The coincidence of Zickendraht coordinates with our collective coordinates depends on 
particle numbers; typically they do not span an invariant manifold. This confirms the well- 
known fact that in general the collective motion in these coordinates does not separate 
rigorously, but only in some adiabatic approximation. 

As we found more general invariant configurations which in turn induce collective coor- 
dinates we may hope that these are useful for approximate considerations for larger particle 
numbers that do include the corresponding invariant manifold in a non-trivial fashion. The 
construction of appropriate coordinates is an open problem. 
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FIGURES 
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FIG. 1. Collective configuration on invariant manifold. Positions are indicated by filled circles 
and momenta by arrows. 






(a) (b) 

FIG. 2. Configurations of eight (a) or six (b) particles in three dimensions that correspond to 
invariant manifolds. Positions are indicated by filled circles. 
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